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This article is the first in a series of three papers investigating the detailed geometry of river net- 
works. Branching networks are a universal structure employed in the distribution and collection 
of material. Large-scale river networks mark an important class of two-dimensional branching 
networks, being not only of intrinsic interest but also a pervasive natural phenomenon. In the 
description of river network structure, scaling laws are uniformly observed. Reported values of 
scaling exponents vary suggesting that no unique set of scaling exponents exists. To improve 
this current understanding of scaling in river networks and to provide a fuller description of 
branching network structure, here we report a theoretical and empirical study of fluctuations 
about and deviations from scaling. We examine data for continent-scale river networks such as 
the Mississippi and the Amazon and draw inspiration from a simple model of directed, random 
networks. We center our investigations on the scaling of the length of sub-basin's dominant 
stream with its area, a characterization of basin shape known as Hack's law. We generalize 
this relationship to a joint probability density and provide observations and explanations of 
deviations from scaling. We show that fluctuations about scaling are substantial and grow with 
system size. We find strong deviations from scaling at small scales which can be explained by 
the existence of linear network structure. At intermediate scales, we find slow drifts in exponent 
values indicating that scaling is only approximately obeyed and that universality remains inde- 
terminate. At large scales, we observe a breakdown in scaling due to decreasing sample space 
and correlations with overall basin shape. The extent of approximate scaling is significantly 
restricted by these deviations and will not be improved by increases in network resolution. 

PACS numbers: 64.60.Ht, 92.40.Fb, 92.40.Gc, 68.70.-Fw 



I. INTRODUCTION 

Networks are intrinsic to a vast number of complex 
forms observed in the natural and man-made world. Net- 
works repeatedly arise in the distribution and sharing of 
information, stresses and materials. Complex networks 
give rise to interesting mathematical and physical proper- 
ties as observed in the Internet , the "small- world" phe- 
nomenon the cardiovascular system ||^, force chains 
in granular media Q] , and the wiring of the brain . 

Branching, hierarchical geometries make up an impor- 
tant subclass of all networks. Our present investigations 
concern the paradigmatic example of river networks. The 
study of river networks, though more general in applica- 
tion, is an integral part of geomorphology, the theory of 
earth surface processes and form. Furthermore, river net- 
works are held to be natural exemplars of allometry, i.e., 
how the dimensions of different parts of a structure scale 
or grow with respect to each other 0, ^, ^, 0. The 
shapes of drainage basins, for example, are reported to 



elongate with increasing basin size ||T^, 13 . 

At present, there is no generally accepted theory 
explaining the origin of this allometric scaling. The fun- 
damental problem is that an equation of motion for ero- 
sion, formulated from first principles, is lacking. The sit- 
uation is somewhat analogous to issues surrounding the 
description of the dynamics of granular media |14, p^. 



noting that erosion is arguably far more complex. Nev- 
ertheless, a number of erosion equations have be en p ro- 
posed ranging from deterministic jl^, [l7|, ^ |l9| to 
stochastic theories ||2^, |2[ |23[^, gj. Each of 
these models attempts to describe how eroding sur- 
faces evolve dynamically. In addition, various heuris- 
tic models of both surface and network evolution also 
exist. Examples include simple lattice-based models of 
erosion |p^, p6| , p7| , p8| , an analogy to invasion per- 
colation p9[, the use of optimality principles and self- 
organized criticality fi , pO|, pT| , and even uncorrelated 
random networks 32|, 33|. Since river networks are 
an essential feature of eroding landscapes, any appropri- 
ate theory of erosoion must yield surfaces with network 
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structures comparable to that of the real world. Howev- 
er, no model of eroding landscapes or even simply of net- 
work evolution unambiguously reproduces the wide range 
of scaling behavior reported for real river networks. 

A considerable problem facing these theories and mod- 
els is that the values of scaling exponents for river net- 
work scaling laws are not precisely known. One of the 
issues we address in this work is universality [p^ p^ . 
Do the scaling exponents of all river networks belong to 
a unique universality class or are there a set of classes 
obtained for various geomorphological conditions? For 
example, theoretical models suggest a variety of expo- 
nent values for networks that are directed versus non- 
directed, created on landscapes with heterogeneous ver- 
sus homogeneous erosivity and so on |3^ . Clearly, 
refined measurements of scaling exponents are imperative 
if we are to be sure of any network belonging to a par- 
ticular universality class. Moreover, given that there is 
no accepted theory derivable from simple physics, more 
detailed phenomenological studies are required. 

Motivated by this situation, we perform here a detailed 
investigation of the scaling properties of river networks. 
We analytically characterize fluctuations about scaling 
showing that they grow with system size. We also report 
significant and ubiquitous deviations from scaling in real 
river networks. This implies surprisingly strong restric- 
tions on the parameter regimes where scaling holds and 
cautions against measurements of exponents that ignore 
such limitations. In the case of the Mississippi basin, 
for example, we find that although our study region 
span four orders of magnitude in length, scaling may be 
deemed valid over no more than 1.5 orders of magnitude. 
Furthermore, we repeatedly find the scaling within these 
bounds to be only approximate and that no exact, single 
exponent can be deduced. We show that scaling breaks 
down at small scales due to the presence of linear basins 
and at large scales due to the inherent discreteness of net- 
work structure and correlations with overall basin shape. 
Significantly, this latter correlation imprints upon river 
network structure the effects and history of geology. 

This paper is the first of a series of three on river- 
network geometry. Having addressed scaling laws in 
the present work, we proceed in second and third arti- 
cles [ p6| , |37| to consider river network structure at a 
more detailed level. In we examine the statistics of 
the "building blocks" of river networks, i.e., segments of 
streams and sub-networks. In particular, we analytically 
connect distributions of various kinds of stream length. 
Part of this material is employed in the present article 
and is a direct generalization of Horton's laws |3^ . In 
the third article , we proceed from the findings of ||3^ 
to characterize how these building blocks fit together. 
Central to this last work is the study of the frequency 
and spatial distributions of tributary branches along the 
length of a stream and is itself a generalization of the 
descriptive picture of Tokunaga @ |l|, |2| . 



II. BASIN ALLOMETRY 

A Hack's law 

In addressing these broader issues of scaling in branch- 
ing networks, we set as our goal to understand the river 
network scaling relationship between basin area a and 
the length I of a basin's main stream: 

l<xa^. (1) 

Known as Hack's law ||ri|, this relation is central to the 
study of scaling in river networks |lj, Hack's expo- 
nent h is empirically found to lie in the range from 0.5 
to 0.7 |ll|, |l|, |l|, @, H, H H, I?! |§. Here, we pos- 
tulate a generalized form of Hack's law that shows good 
agreement with data from real world networks. 

We focus on Hack's law because of its intrinsic interest 
and also because many interrelationships between a large 
number of scaling laws are known and only a small subset 
are understood to be independent ]T^ , p9| . Thus, our 
results for Hack's law will be in principle extendable to 
other scaling laws. With this in mind, we will also discuss 
probability densities of stream length and drainage area. 

Hack's law is stated rather loosely in equation (|l|) and 
implicitly involves some type of averaging which needs 
to be made exphcit. It is most usually considered to be 
the relationship between mean main stream length and 
drainage area, i.e., 

(0 cx a*". (2) 

Here, (•) denotes ensemble average and {I) = {l{a)) is the 
mean main stream length of all basins of area a. Typi- 
cally, one performs regression analysis on log (l) against 
logo to obtain the exponent h. 

B Fluctuations and deviations 

In seeking to understand Hack's law, we are natural- 
ly led to wonder about the underlying distribution that 
gives rise to this mean relationship. By considering fluc- 
tuations, we begin to see Hack's law as an expression of 
basin morphology. What shapes of basins characterized 
by (a, I) are possible and with what probability do they 
occur? 

An important point here is that Hack's law does not 
exactly specify basin shapes. An additional connection to 
Euclidean dimensions of the basin is required. We may 
think of a basin's longitudinal length and its width 
L±. The main stream length I is reported to scale with 
L|| as 

/ cx Ll (3) 

where typically 1.0 < d < 1.1, [|l2|, |5^. Hence, we have 
a oc l^^'^ oc Ly^''. All other relevant scaling laws expo- 
nents can be related to the pair of exponents (d, h) which 
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FIG. 1: The Full Hack distribution for the Kansas, (a), 
and Mississippi, (b), river basins. For each value of a, the 
distribution has been normalized along the I direction by 
max, Pia,l) |§. 



therefore characterize the universality class of a river net- 
work |]l^, U d/h — 2 we have that basins are self- 
similar whereas ii d/h < 2, we have that basins are elon- 
gating. So, while Hack's law gives a sense of basin allom- 
etry, the fractal properties of main stream lengths need 



also be known in order to properly quantify the scaling 
of basin shape. 

In addition to fluctuations, complementary insights are 
provided by the observation and understanding of devi- 
ations from scaling. We are thus interested in discerning 
the regularities and quirks of the joint probability dis- 
tribution P{a,l). We will refer to P{a,l) as the Hack 
distribution. 

Hack distributions for the Kansas river basin and 
the Mississippi river basin are given in Figures |^(a) 
and|^(b). Fluctuations about and deviations from scaling 
are immediately evident for the Kansas and to a lesser 
extent for the Mississippi. The first section of the paper 
will propose and derive analytic forms for the Hack distri- 
bution under the assumption of uniform scaling with no 
deviations. Here, as well as in the following two papers 
of this series ^ , we will motivate our results with a 
random network model originally due to Scheidegger [ p3[ . 

We then expand our discussion to consider deviations 
from exact scaling. In the case of the Kansas river, a 
striking example of deviations from scaling is the linear 
branch separated from the body of the main distribu- 
tion shown in Figure [^(a) . This feature is less prominent 
in the lower resolution Mississippi data. Note that this 
linear branch is not an artifact of the measurement tech- 
nique or data set used. This will be explained in our dis- 
cussion of deviations at small scales in the paper's second 
section. 

We then consider the more subtle deviations associated 
with intermediate scales. At first inspection, the scaling 
appears to be robust. However, we find gradual drifts in 
"exponents" that prevent us from identifying a precise 
value of h and hence a corresponding universality class. 

Both distributions also show breakdowns in scaling for 
large areas and stream lengths and this is addressed in 
the final part of our section on deviations. The reason for 
such deviations is partly due to the decrease in number 
of samples and hence self-averaging, as area and stream 
lengths are increased. However, we will show that the 
direction of the deviations depends on the overall basin 
shape. We will quantify the extent to which such devi- 
ations can occur and the effect that they have on mea- 
surements of Hack's exponent h. 

Throughout the paper, we will return to the Hack dis- 
tributions for the Kansas and the Mississippi rivers as 
well as data obtained for the Amazon, the Nile and the 
Congo rivers. 

III. FLUCTUATIONS: AN ANALYTIC FORM 
FOR THE HACK DISTRIBUTION 

To provide some insight into the nature of the under- 
lying Hack distribution, we present a line of reasoning 
that will build up from Hack's law to a scaling form of 
P{a,l). First let us assume for the present discussion of 
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fluctuations that an exact form of Hack's law holds: 

(/) = 6la'' (4) 

where we have introduced the coefficient 9 which we dis- 
cuss fully later on. Now, since Hack's law is a power law, 
it is reasonable to postulate a generalization of the form 



P{l\a) 



1 



I 



(5) 



The prefactor l/a'* provides the correct normalization 
and Fi is the "scaling function" we hope to understand. 
The above will be our notation for all conditional prob- 
abilities. Implicit in equation is the assumption that 
all moments and the distribution itself also scale. For 
example, the qth moment of P{1\ a) is 



which implies 



(Z«(a)) cx a"''. 



(6) 



(7) 



where S K. Also, for the distribution P{1 \ a) it follows 
from equation (^) that 



k''P{lk^ I ak) 



-F, 



Ik'' 



Pil\a). (8) 



We note that previous investigations of Hack's law jlj, 
consider the generalization in equation (^) . Rigon et 
al. also examine the behavior of the moments of the 
distribution P(l \ a) for real networks. Here, we will go 
further to characterize the full distribution P(a, I) as well 
as both P{1 1 a) and P(a | /). Along these lines, Rigon et 
al. suggest that the function Fi{x) is a "finite-size" 
scaling function analogous to those found in statistical 
mechanics, i.e.: Fi{x) — > as a; ^ oo and Fi{x) — > c as 
a; — > 0. However, as we will detail below, the restrictions 
on Fiix) can be made stronger and we will postulate a 
simple Gaussian form. More generally, Fiix) should be 
a unimodal distribution that is non-zero for an interval 
[a;i,a;2] where xi > 0. This is so because for any given 
fixed basin area a, there is a minimum and maximum I 
beyond which no basin exists. This is also clear upon 
inspection of Figures |l|(a) and ^(b). 

We observe that neither drainage area nor main stream 
length possess any obvious features so as to be deemed 
the independent variable. Hence, we can also view Hack's 
law as its inversion (a) oc Note that the constant of 
proportionality is not necessarily 0^^'' and is dependent 
on the nature of the full Hack distribution. We thus have 
another scaling ansatz as per equation ^ 



P{a\l) = l/f/''Faia/l'/''). 



(9) 



The conditional probabilities P{1 \ a) and P{a \ I) are 
related to the joint probability distribution as 



where P(l) and P{a) are the probability densities of main 
stream length and area. These distributions are in turn 
observed to be power laws both in real world networks 
and models ^ ^ : 



P{a) - Naa 



and 



P(0 ~ Nil-\ 



(11) 



where Na and Ni are appropriate prefactors and the tilde 
indicates asymptotic agreement between both sides for 
large values of the argument. Furthermore, the expo- 
nents T and 7 are related to Hack's exponent h via the 
scaling relations p9] 



2-h 



and 



7 = 1/h. 



(12) 



Equations (|), (|), (|ol), (|rT|), and (H) combine to give 
us two forms for P{a, I), 

P{a, I) = \/a^F{l/a'') = l/l'^f'^Gia/l^'''), (13) 

where x^'^ F{x^'') — G(x) and, equivalently, F{y) — 



IV. RANDOM DIRECTED NETWORKS 

We will use results from the Scheidegger model to 
motivate the forms of these distributions. In doing so, 
we will also connect with some problems in the theory of 
random walks. 



A 
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P{a, I) = P{a)P{l I a) = P{l)P{a \ I), 



(10) 



FIG. 2: Scheidegger's model of random, directed net- 
works. Flow is down the page and at each site, stream 
flow is randomly chosen to be in one of the two downward 
diagonals. Stream paths and basin boundaries are thus 
discrete random walks. 



Scheidegger's model of river networks is defined on a 
triangular lattice as indicated by Figure |[ Flow in the 
figure is directed down the page. At each site, the stream 
fiow direction is randomly chosen between the two diag- 
onal directions shown. Periodic boundary conditions are 
applied in all of our simulations. Each site locally drains 
an area of a^, where the lattice unit a is the distance 
between neighboring sites, and each segment of stream 
has a length a. For simplicity, we will take a to be unity. 
We note that connections exist between the Scheideg- 
ger model and models of particle aggregation pU p3[. 



5 



Abelian sandpiles |5^ and limiting cases of force 

chain models in granular media 

Since Scheidegger's model is based on random flow 
directions, the Hack distributions have simple interpre- 
tations. The boundaries of drainage basins in the model 
are random walks. Understanding Hack's law therefore 
amounts to understanding the first collision time of two 
random walks that share the same origin in one dimen- 
sion. If we subtract the graph of one walk from the other, 
we see that the latter problem is itself equivalent to the 
first return problem of random walks [ p7| . 

Many facets of the first return problem are well under- 
stood. In particular, the probability of n, the number of 
steps taken by a random walk until it first returns to the 
origin, is asymptotically given by 

Pin) ^ -^n-^'\ (14) 

But this number of steps is also the length of the basin 
I. Therefore, we have 



P{1) 



(15) 



where because we are considering the difference of two 
walks, we use P{1) = P(n/2)|„=;. Also, we have found 
the pref actor Ni = 2/ti. 

We thus have that 7 = 3/2 for the Scheidegger model. 
The scaling relations of equation ( p^ ) then give /i = 2/3 
and T = 4/3. The value of h is also readily obtained 
by noting that the typical area of a basin of length I is 
a oc I - V-/"^ = P/'^ = V-/^ since the boundaries are random 
walks. 



P(a/ ^'^) 




FIG. 3: Cross-Sectional scaling functions of the Hack 
distribution for the Scheidegger model with lattice con- 
stant equal to unity. Both distributions are normalized. 
The right distribution is for I fixed and a varying and 
is postulated to be a normal distribution. The left dis- 
tribution for a fixed and I varying and is a form of an 
inverse Gaussian Q. The data used was obtained for 
all sites with with I > 100 and a > 500 respectively. 
Each distribution was obtained from ten realizations of 
the Scheidegger model on a 10'* x 10* lattice. 



V. AREA-LENGTH DISTRIBUTION FOR 
RANDOM, DIRECTED NETWORKS 

Something that is less well studied is the joint distri- 
bution of the area enclosed by a random walk and the 
number of steps to its first return. In terms of the Schei- 
degger model, this is precisely the Hack distribution. 

We motivate some general results based on observa- 
tions of the Scheidegger model. Figure || shows the nor- 
malized distributions P(a^^"^/^) and P{la^^/^) as derived 
from simulations of the model. Given the scaling ansatzes 
for P{1 1 a) and P{a \ I) in equations (^) and (||), we see 
that P{y = la-^/'^) = Fi{y) and P{x = al-^^^) = Fa(x). 

Note that we have already used Hack's law for the 
Scheidegger model with /i = 2/3 to obtain these distri- 
butions. The results are for ten realizations of the model 
on a lO'' by lO'' lattice, taking 10^ samples from each of 
the ten instances. For P{a \ I), only sites where I > 100 
were taken, and similarly, for P{l\a), only sites where 
a > 500 were included in the histogram. 

We postulate that the distribution P{y — la~^^^) is a 
Gaussian having the form 

P{1 1 a) = ^=^— exp{-(?a-2/3 _ ef /2r^'} (16) 



We estimate the mean of Fi to be 6* ~ 1.675 (this is 
the same 6 as found in equation (^)) and the standard 
deviation to be 77 ~ 0.321. The fit is shown in Figure ^ 
as a solid line. The above equation agrees with the form 
of the scaling ansatz of equation (|^) and we now have 
the assertion that Fi is a Gaussian defined by the two 
parameters 9 and 77. 

Note that the 9 and 77 are coefficients for the actual 
mean and standard deviation. In other words, for fixed 
a, the mean of P{1 \ a) is 90^^^ and its standard deviation 
is Tjo?/^. Having observed their context, we will refer to 
9 and 77 as the Hack mean coefficient and Hack standard 
deviation coefficient. 

From this starting point we can create P{a, I) and 
P{a I Z), the latter providing a useful test. Since P{a) ^ 
Naa^'^ — Naa~'^/'^, as per equation (O), we have 



P(a,0 



a4/3 V2^a2/3,7 
N, 



exp{-(/a-2/3_^)2/2^2|^ 
exp{-(/a-2/3_ 0)2/2772}. (17) 



27ra2/3?7 



As expected, we observe the form of equation (^7|) to be 
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in accordance with that of equation ( p^ . Note that the 
scahng function F (and equivalently G) is defined by the 
three parameters 0, rj and Na, the latter of which may 
be determined in terms of the former as we will show 
below. Also, since we expect all scaling functions to be 
only asymptotically correct, we cannot use equation (|l7| ) 
to find an expression for the normalization Na- Equa- 
tion ( ITtI) ceases to be valid for small a and I. However, 
we will be able to do so once we have P{a \ I) since we 
are able to presume I is large and therefore that the scal- 
ing form is exact. Using equation (|l^) and the fact that 
P{a 1 1) = P(a, I)/ P{1) from equation ( p^ we then have 



P[a\l) 



,3/2 



exp{~G/a'/'-^)V27?'}, 



^|^exp{-a/a^/^-^)V2.n, 



1 

13/2 23/2^ 

xexp{-((a/;3/2)-2/3_0)2/2^2|_ 



(18) 



In rearranging the expression of P{a \ I), we have made 
clear that its form matches that of equation (||). 

A closed form expression for the normalization factor 
Na may now be determined by employing the fact that 
JZodaP{a\l)^l. 



1 = I daP{a\l), 

la=Q 

da N 



{a/l 



3/2N 



,^0 ;3/2 23/2^ 

xexp{-((a/;3/2)-2/3_0)2/2^2|^ 

Na^ 



25/2 



Auu^'^eyiY>{-{u-ef/27f}, (19) 



where we have used the substitution ajl^l"^ = u ^/2 and 
hence also Z~^/2 = (— 3/2)-u~^/2 \Ye therefore have 



Na^ 



25/2, 



/3n 



(20) 

We may thus write down all of the scaling functions 
Fi, Fa, F and G for the Scheidegger model: 



Fiiz) 

Fa{z) 

F{z) 
G{z) 



1 



27r?7 



^" exp{-(z- 0)2/2772}, and 



2711] 
2Nr 



(21) 
(22) 
(23) 



0F23/2,7 



z-2exp{-(z-2/3 _6))2/2,72|_ (24) 



Recall that all of these forms rest on the assumption 
that Fi{z) is a Gaussian. In order to check this assump- 
tion, we return to Figure |^. The empirical distribution 



P{z = a/P/"^) is shown on the left marked with cir- 
cles. The solid line through these points is Fa{z) as given 
above in equation ( pT| ) . There is an excellent match so we 
may be confident about our proposed form for Fa{z). We 
note that the function Fa{z) may be thought of as a frac- 
tional inverse Gaussian distribution, the inverse Gaussian 
being a well known distribution arisin g in the study of 
first passage times for random walks [p7| . It is worth 
contemplating the peculiar form of P{a \ I) in terms of 
first return random walks. Here, we have been able to 
postulate the functional form of the distribution of areas 
bound by random walks that first return after n steps. 
If one could understand the origin of the Gaussian and 
find analytic expressions for 9 and 77, then the problem 
would be fully solved. 



VI. AREA-LENGTH DISTRIBUTION 
EXTENDED TO REAL NETWORKS 

We now seek to extend these results for Scheidegger's 
model to real world networks. We will look for the same 
functional forms for the Hack distributions that we have 
found above. The conditional probability distributions 
pertaining to Hack's law take on the forms 



P{l\a) = l/a-''Fi{la-'') 



and 

P{a\l)^l-^"'Fa{al-^"') 



2TTr] 



exp{-(Za-"-(?)V2ry2}, 
(25) 



1-1/ h 



Na 



2tt Ni-q 



(a/-i/'')-2 exp{-((ari/'')-'' - ef/2jf}, 



(26) 



and the full Hack distribution is given by 

P{a,l) = a^Fila-'^) = l-'^'''G{al-^/'') 

Na 



2'Kria^ 



exp{-(/a-"- 0)2/2772}. (27) 



with Na determined by equation (pO|). The three param- 
eters /i, the Hack exponent, 9, the Hack mean coeffi- 
cient, and 77, the Hack standard deviation coefRcient, are 
in principle landscape dependent. Furthermore, in 77 we 
have a basic measure of fluctuations in the morphology 
of basins. 

Figures ^ and || present Hack scaling functions for the 
Mississippi and Nile river basins. These Figures is to be 
compared with the results for the Scheidegger model in 
Figure |3[ 

For both rivers, Hack's exponent h was determined 
first from a stream ordering analysis (we discuss stream 
ordering later in Section ^ . Estimates of the parameters 
9 and 77 were then made using the scaling function Fi 
presuming a Gaussian form. 

We observe the Gaussian fit for the Mississippi is more 
satisfactory than that for the Nile. These fits are not 
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12 3 4 

FIG. 4: Cross-sectional scaling functions of the Hack 
distribution for the Mississippi. The estimates used for 
Hack's exponent are h — 0.55 and h = 0.50 respecti vely , 
the determination of which is discussed in section [VI]| . 
The fits indicated by the smooth curves to the data 
are made as per the Scheidegger model in Figure |^ and 
according to equations (^) and (^). The values of the 
Hack mean coefficient and standard deviation coefficient 
are estimated to be 6* ~ 0.80 and rj ~ 0.20. 



rigorously made because even though we have chosen 
data ranges where deviations (which we address in the 
following section) are minimal, deviations from scaling 
do still skew the distributions. The specific ranges used 
to obtain Fi and Fa respectively are for the Mississippi: 
8.5 < logio ^ < 9.5 and 4.75 < log^Q I < 6, and for the 
Nile: 9 < log^Q a < 11 and 5.5 < log^o I < 6 (areas are in 
km^ and lengths km). Furthermore, we observe that the 
estimate of h has an effect in the resulting forms of Fi 
and Fa- Nevertheless, here we are attempting to capture 
the essence of the generalized form of Hack's law in real 
networks. 

We then use the parameters h, 9 and rj and equa- 
tion ( |26| ) to construct our theoretical Fa, the smooth 
curves in Figures || and ||. As for the Scheidegger model 
data in Figure H, we see in both examples approximate 
agreement between the measured Fa and the one predict- 
ed from the form of F/. Table | shows estimates of h, 9 
and r] for the five major river basins studied. 

Given our reservations about the precision of these val- 
ues of 9 and 77, we are nevertheless able to make qualita- 
tive distinctions. Recalling that (/) = 9a^, we see that, 
for fixed h, higher values of 9 indicate relatively longer 
stream lengths for a given area and hence longer and thin- 




FIG. 5: Cross-sectional scaling functions of the Hack 
distribution for the Nile. The top and right axes corre- 
spond to Fa and the bottom and left to Fi. The Hack 
exponent used is h = 0.50 and the values of the Hack 
mean coefficient and standard deviation coefficient are 
estimate to be 6* ~ 2.45 and 77 ~ 0.50 
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Mississippi 


0.80 


0.20 


0.55 


Amazon 


1.90 


0.35 


0.52 


Nile 


2.45 


0.50 


0.50 


Kansas 


0.70 


0.15 


0.57 


Congo 


0.89 


0.18 


0.54 



TABLE I: Estimates of Hack distribution parameters 
for real river networks. The scaling exponent h is Hack's 
exponent. The parameters 9 and 77 are coefficients of the 
mean and standard deviation of the conditional proba- 
bility density function P{1 \ a) and are fully discussed in 
the text @. 



ner basins. The results therefore suggest the Nile, and 
to a lesser degree the Amazon, have basins with thinner 
profiles than the Congo and, in particular, the Mississippi 
and Kansas. This seems not unreasonable since the Nile 
is a strongly directed network constrained within a rel- 
atively narrow overall shape. This is somewhat in spite 
of the fact that the shape of an overall river basin is not 
necessarily related to its internal basin morphology, an 
observation we will address later in Section 

The importance of 9 is tempered by the value of h. 
Hack's exponent affects not only the absolute measure of 
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stream length for a given area but also how basin shapes 
change with increasing area. So, in the case of the Kansas 
the higher value of h suggests basin profiles thin with 
increasing size. This is in keeping with overall direct- 
edness of the network. Note that our measurements of 
the fractal dimension d of stream lengths for the Kansas 
place it to be d = 1.04 ± 0.02. Therefore, d/h ~ 1.9 < 2 
and elongation is still expected when we factor in the 
scaling of I with L||. 

The Nile and Amazon also all have relatively high 77 
indicating greater fluctuations in basin shape. In compar- 
ison, the Mississippi, Kansas and Congo appear to have 
less variation. Note that the variability of the Kansas is 
in reasonable agreement with that of the whole Missis- 
sippi river network for which it is a sub-basin. 

Finally, regardless of the actual form of the distribution 
underlying Hack's law, fluctuations are always present 
and an estimate of their extent is an important measure- 
ment. Thus, the Hack mean and standard deviation coef- 
ficients, 9 and 77, are suggested to be of sufficient worth 
so as to be included with any measurement of the Hack 
exponent h. 



VII. DEVIATIONS FROM SCALING 

In generalizing Hack's law, we have sought out regions 
of robust scaling, discarding ranges where deviations 
become prominent. We now bring our attention to the 
nature of the deviations themselves. 

We observe three major classes of deviations which we 
will define by the scales at which they occur: small, inter- 
mediate and large. Throughout the following sections we 
primarily consider deviations from the mean version of 
Hack's law, (l) = Oa^, given in equation (^. Much of 
the understanding we gain from this will be extendable 
to deviations for higher moments. 

To provide an overview of what follows, examples of 
mean-Hack distributions for the Kansas river and the 
Mississippi are shown in Figures ||(a) and ^(b). Hack's 
law for the Kansas river exhibits a marked deviation 
for small areas, starting with a near linear relationship 
between stream length and area. A long crossover region 
of several orders of magnitude in area then leads to an 
intermediate scaling regime wherein we attempt to deter- 
mine the Hack exponent h. In doing so, we show that 
such regions of robust scaling are surprisingly limited 
for river network quantities. Moreover, we observe that, 
where present, scaling is only approximate and that no 
exact exponents can be ascribed to the networks we study 
here. It follows that the identification of universality 
classes based on empirical evidence is a hazardous step. 

Finally, the approximate scaling of this intermediate 
region then gives way to a break down in scaling at larger 
scales due to low sampling and correlations with basin 
shape. The same deviations are present in the relatively 
coarse-grained Mississippi data but are less pronounced. 



(a) 




4 5 6 7 8 9 10 11 12 




8 9 10 11 12 13 



FIG. 6: The mean version of Hack's law for the Kansas, 
(a), and the Mississippi, (b). The units of lengths and 
areas are meters and square meters. These are calculat- 
ed from the full Hack distributions shown in Figures |^(a) 
and 0(b) by finding (/) for each value of basin area a. 
Area samples are taken every 0.02 orders of magnitude 
in logarithmic space. The upper dashed lines represent a 
slope of unity in both plots and the lower lines the Hack 
exponents 0.57 and 0.55 for the Kansas and Mississippi 
respectively. For the Kansas, there is a clear deviation 
for small area which rolls over into a region of very slow- 
ly changing derivative before breaking up at large scales. 
Deviations from scaling are present for the lower resolu- 
tion dataset of the Mississippi but to a lesser extent. 



9 




3.5 4 4.5 5 5.5 6 

FIG. 7: The linearity of Hack's law at small scales for 
the Kansas river. The linear regime enters a crossover 
region after almost 1.5 orders of magnitude in area. This 
is an expanded detail of the mean Hack's law given in 
Figure ^(a) and areas and lengths are in square meters 
and meters. 



VIII. DEVIATIONS AT SMALL SCALES 

At small scales, we find the mean-Hack distribution to 
follow a linear relationship, i.e., I (x a^. This feature is 
most evident for the Kansas river as shown in Figure |^(a) 
and in more detail in Figure 0. The linear regime per- 
sists for nearly 1.5 orders of magnitude in basin area. To 
a lesser extent, the same trend is apparent in the Missis- 
sippi data. Figure ^(b). 

Returning to the full Hack distribution of Figures |^(a) 
and 0(b), we begin to see the origin of this linear regime. 
In both instances, a linear branch separates from the 
body of the main distribution. Since / cannot grow faster 
than , the linear branch marks an upper bound on the 
extent of the distribution in {a, I) coordinates. When 
averaged to give the mean-Hack distribution, this linear 
data dominates the result for small scales. 

We find this branch evident in all Hack distributions. 
It is not an artifact of resolution and in fact becomes 
more pronounced with increased map precision. The ori- 
gin of this linear branch is simple: data points along the 
branch correspond to positions in narrow sub-networks, 
i.e., long, thin "valleys." To understand the separation 
of this linear branch from the main body of the distribu- 
tion, consider Figure ^ which depicts a stream draining 
such a valley with length and area (Zi,ai) that meets a 




ai a2 as 



FIG. 8: Origin of the linear branch in the Hack distri- 
bution. The sub-basins labeled 1 and 2 depicted on the 
left have areas ai and 02 and lengths h and ^2- These 
sub-basins combine to form a basin of area — ai + 02 
and main stream length ^3 = maxZi,/2 = h- Since sub- 
basin 1 is a linear sub- network (a valley) the pair (ai, h) 
lie along the linear branch of the Hack distribution as 
shown on the right. Points along the main stream of 
sub-basin 1 lie along the dashed line leading to the coor- 
dinate (ai, ^i). On combining with the second basin, the 
jump in the resultant area creates a jump from (oi, h to 
(as, I2) in the main body of the Hack distribution. 



stream from a basin with characteristics (^2,ci2)- The 
area and length of the basin formed at this junction 
is thus (a3,^3) — (ai -|- a2, max(Zi, Z2)) (in the Figure, 
max(Zi,Z2) = h)- The greater jump in area moves the 
point across into the main body of the Hack distribu- 
tion, creating the separation of the linear branch. 

In fact, the full Hack distribution is itself comprised 
of many such linear segments. As in the above example. 
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until a stream does not meet any streams of comparable 
size, then its area and length will roughly increase in 
linear fashion. When it does meet such a stream, there is 
a jump in area and the trace of a new linear segment is 
started in the distribution. We will see this most clearly 
later on when we study deviations at large scales. 

For very fine scale maps, on the order of meters, we 
might expect to pick up the scale of the unchannelized, 
convex regions of a landscape, i.e., "hillslopes" |Q. This 
length scale represents the typical separation of branches 
at a network's finest scale. The computation of stream 
networks for these hillslope regions would result in large- 
ly non-convergent (divergent or parallel) flow. There- 
fore, we would have linear "basins" that would in theory 
contribute to the linear branch we observe |0. Poten- 
tially, the crossover in Hack's law could be used as a 
determinant of hillslope scale, a crucial parameter in geo- 
morphology However, when long, thin network 

structures are present in a network, this hillslope scale is 
masked by their contribution. 

Whether because of the hillslope scale or linear net- 
work structure, we see that at small scales. Hack's law 
will show a crossover in scaling from /i = 1 to a low- 
er exponent. The crossover's position depends on the 
extent of linear basins in the network. For example, 
in the Kansas River basin, the crossover occurs when 
{I, a) = (4 X lO'^m, lO^m^). Since increased map resolu- 
tion can only increase measures of length, the crossover's 
position must occur at least at such a length scale which 
may be many orders of magnitude greater than the scale 
of the map. 

However, the measurement of the area of such lin- 
ear basins will potentially grow with coarse-graining. 
Note that for the Mississippi mean-Hack distribution, 
the crossover begins around a = 10^-^ m^ whereas for 
the Kansas, the crossover initiates near a — 10^'^ m^ but 
the ends of the crossovers in both cases appear to agree, 
occurring at around a — 10^ m^. Continued coarse- 
graining will of course eventually destroy all statistics 
and introduce spurious deviations. Nevertheless, we see 
here that the deviation which would only be suggested in 
the Mississippi data is well confirmed in the finer-grain 
Kansas data. 



IX. DEVIATIONS AT INTERMEDIATE SCALES 

As basin area increases, we move out of the linear 
regime, observing a crossover to what would be consid- 
ered the normal scaling region of Hack's law. We detail 
our attempts to measure the Hack exponents for the 
Kansas and Mississippi examples. Rather than relying 
solely on a single regression on a mean-Hack distribu- 
tion, we employ a more precise technique that examines 
the distribution's derivative. As we will show, we will not 
be able to find a definite value for the Hack exponent in 
either case, an important result in our efforts to deter- 



mine whether or not river networks belong to specific 
universality classes. 

To determine h, we consider Hack's law (equation 
explicitly in logarithmic coordinates, 



logio (Ka)> = logio ^ + logio a- 



(28) 



The derivative of this equation with respect to logj^g ^ 
then gives Hack's exponent as a function of area. 



h{a) 



d logio a 



logio (Ka))- 



(29) 



We may think of h{a) as a "local Hack exponent." Note 
that non-constant trends in h{a) indicate scaling does not 
hold. We calculate the discrete derivative as above for the 
Kansas and Mississippi. We smooth the data by taking 
running averages with varying window sizes of n sam- 
ples, the results for n = 32 being shown in Figures ||(a) 
and ^(b) where the spacing of log a is 0.02 orders of mag- 
nitude. Thus, the running averages for the figures are 
taken over corresponding area ranges 0.64 orders of mag- 
nitude. 

Now, if the scaling law in question is truly a scaling 
law, the above type of derivative will fluctuate around a 
constant value of exponent over several orders of mag- 
nitude. With increasing n, these fluctuations will neces- 
sarily decrease and we should see the derivative holding 
steady around the exponent's value. 

At first glance, we notice considerable variation in h{a) 
for both data sets with the Kansas standing out. Fluctu- 
ations are reduced with increasing n but we observe con- 
tinuous variation of the local Hack exponent with area. 
For the example of the Kansas, the linear regime and 
ensuing crossover appear as a steep rise followed by a 
drop and then another rise during all of which the local 
Hack exponent moves well below 1/2. 

It is after these small scale fiuctuations that we would 
expect to find Hack's exponent. For the Kansas river 
data, we see the derivative gradually climbs for all values 
of n before reaching the end of the intermediate regime 
where the putative scaling breaks down altogether. 

For the Kansas show in Figure ||(a), the dashed fine 
represents h = 0.57, our estimate of Hack's exponent 
from simple regression on the mean-Hack distribution of 
Figure ^(a). For the regression calculation, the interme- 
diate region was identified from the figure to be 10^ < 
a < 10^° m^. We see from the smoothed derivative in 
Figure |9|(a) that the value h = 0.57 is not precise. After 
the crossover from the linear region has been complet- 
ed, we observe a slow rise from h ~ 0.54 to ft, ~ 0.63. 
Thus, the local Hack exponent h{a) gradually climbs 
above h = 0.57 rather than fluctuate around it. 

A similar slow change in h(a) is observed for the Missis- 
sippi data. We see in Figure ||(b) a gradual rise and then 
fall in h{a). The dashed line here represents h = 0.55, 
the value of which was determined from Figure ^(b) using 
regression on the range 10® < a < 10^^'^ m^. The range 
of h{a) is roughly [0.52,0.58]. Again, while h = 0.55 
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FIG. 9: Variation in Hack's exponent for the Kansas, 
(a), and the Mississippi, (b). The area a is in m^. The 
plots are derivatives of the mean-Hack distributions given 
in Figures ^(a) and ^(b). Both derivatives have been 
smoothed by taking running averages over 0.64 orders 
of magnitude in a. For the Kansas, the dashed hue is 
set dX h = 0.57, Hack's exponent estimated via simple 
regression analysis for points with 10'' < a < 10^° m^. 
The local exponent is seen to gradually rise through the 
h = 0.57 level indicating scaling is not robust. For the 
Mississippi in (b), the dashed line is a Hack exponent of 
0.55 calculated from regression on data in the interval 
10^ < a < 10^^'^ m^. The local Hack exponent is seen to 
gradually rise and fall about this value. 



approximates the derivative throughout this intermedi- 
ate range of Hack's law, we cannot claim it to be a precise 
value. 

We observe the same drifts in h{a) in other datasets 
and for varying window size n of the running aver- 
age. The results suggest that we cannot assign specific 
Hack exponents to these river networks and are therefore 
unable to even consider what might be an appropriate 
universality class. The value of h obtained by regression 
analysis is clearly sensitive to the the range of a used. 
Furthermore, these results indicate that we should main- 
tain healthy reservations about the exact values of other 
reported exponents. 



X. DEVIATIONS AT LARGE SCALES 

We turn now to deviations from Hack's law at large 
scales. As we move beyond the intermediate region of 
approximate scaling, fluctuations in h{a) begin to grow 
rapidly. This is clear on inspection of the derivatives 
of Hack's law in Figures ||(a) and |9|(b). There are two 
main factors conspiring to drive these fluctuations up. 
The first is that the number of samples of sub-basins 
with area a decays algebraically in a. This is just 
the observation that P{a) oc as per equation (|ll]). 
The second factor is that fluctuations in I and a are 
on the order of the parameters themselves. This fol- 
lows from our generalization of Hack's law which shows, 
for example, that the moments (^') of P{1 \ a) grow like 



qh 



Thus, the standard deviation grows like the mean: 



cf{l) = {{P) - (Z)2)i/2 



A Stream ordering and Horton's laws 

So as to understand these large scale deviations from 
Hack's law, we need to examine network structure in 
depth. One way to do this is by using Horton-Strahler 
stream ordering 1 38 , 63] a nd a generalization of the well- 
known Horton's laws [pSl, |6^ , |65| , ^ ^ . This will natu- 
rally allow us to deal with the discrete nature of a network 
that is most apparent at large scales. 

Stream ordering discretizes a network into a set of 
stream segments (or, equivalently, a set of nested basins) 
by an iterative pruning. Source streams (i.e., those with- 
out tributaries) are designated as stream segments of 
order uj — \. These are removed from the network and 
the new source streams are then labelled as order oj = 2 
stream segments. The process is repeated until a single 
stream segment of order w = is left and the basin itself 
is defined to be of order Vl. 

Natural metrics for an ordered river network are rii^, 
the number of order lo stream segments (or basins), a^;, 
the average area of order oj basins, l^^, the average main 



is] 

stream length of order uj basins, and l^s , the aver- 
age length of order uj stream segments. Horton's laws 
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FIG. 10: Hack distribution for the Mississippi plotted in 
linear space. The area units a is 10^^ and length I is 
10^ m. The discreteness of the basin structure is clearly 
indicated by the isolated, linear fragments. The inset is 
a blow-up of the box on the main graph. 



state that these quantities change regularly from order 
to order, i.e.. 



Rn and 



X. 



UJ+l 



Rx, 



(30) 



where X = a, I or Note that all ratios are defined to 
be greater than unity since areas and lengths increase but 
number decreases. Also, there are only two independent 



ratios since Ra = Rn and Ri = R 



Horton's laws 



mean that stream-order quantities change exponentially 
with order. For example, (pO|) gives that cx (i?/)". 



B Discrete version of Hack's law 

Returning to Hack's law, we examine its large scale 
fluctuations with the help of stream ordering. We are 
interested in the size of these fluctuations and also how 
they might correlate with the overall shape of a basin. 
First, we note that the structure of the network at large 
scales is explicitly discrete. Figure |l^ demonstrates this 
by plotting the distribution of (a, I) without the usual 
logarithmic transformation. Hack's law is seen to be 
composed of linear fragments. As explained above in 
Figure ^, areas and length increase in proportion to each 
other along streams where no major tributaries enter. As 



soon as a stream does combine with a comparable one, a 
jump in drainage area occurs. Thus, we see in Figure |lO| 
isolated linear segments which upon ending at a point 
(oi, h) begin again at (oi -I- a2, ^i), i.e., the main stream 
length stays the same but the area is shifted. 

We consider a stream ordering version of Hack's law 
given by the points {au),lu})- The scaling of these data 
points is equivalent to scaling in the usual Hack's law. 
Also, given Horton's laws, it follows that h = Ini?;/ lni?„ 
(using Ra = -Rn). Along the lines of the derivative we 
introduced to study intermediate scale fluctuations in 
equation (E9h, we have here an order-based difference: 



log lui/lul 

log a^/K 



(31) 



We can further extend this definition to differences 
between non-adjacent orders: 



log L/L 

log d^/a^ 



(32) 



This type of difference, where to' < may be best 
thought of as a measure of trends rather than an approx- 
imate discrete derivative. 

Using these discrete differences, we examine two fea- 
tures of the order-based versions of Hack's law. First we 
consider correlations between large scale deviations with- 
in an individual basin and second, correlations between 
overall deviations and basin shape. For the latter, we 
will also consider deviations as they move back into the 
intermediate scale. This will help to explain the gradual 
deviations from scaling we have observed at intermediate 
scales. 

Since deviations at large scales are reflective of only a 
few basins, we require an ensemble of basins to provide 
sufficient statistics. As an example of such an ensemble, 
we take the set of order $7 = 7 basins of the Mississippi 
basin. For the dataset used here where the overall basin 
itself is of order = 11, we have 104 order 17 = 7 sub- 
basins. The Horton averages for these basins are 07 ~ 
16600 km^, It ~ 350 km, and L7 ~ 210 km. 

For each basin, we first calculate the Horton averages 
{dujjuj)- We then compute h^^^^i, the Hack difference 
given in equation ( |3l|) . To give a rough picture of what 
is observed, Figurc]ll| shows a scatter plot of h^^^-i for 
all order = 7 basins. Note the increase in fluctua- 
tions with increasing lj. This increase is qualitatively 
consistent with the smooth versions found in the single 
basin examples of Figures ^(a) and ||(b). In part, less 
self-averaging for larger ld results in a greater spread in 
this discrete derivative. However, as we will show, these 
fluctuations are also correlated with fluctuations in basin 
shape. 
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FIG. 11: Differences of the stream order-based version of 
Hack's law for 104 order fl = 7 basins of the Mississippi 
(compare the continous versions given in Figure |^) The 
plots are overlaid to give a sense of the increase in fluctu- 
ations of the local Hack exponent h^.uj-i with increasing 
order to. The clusters correspond to w = 1, 2, . . . , 7, mov- 
ing across from left to right. 



FIG. 12: A comparison of the stream order Hack deriva- 
tives hn,n-i and /ia-1,0-2 for each of the 104 order ft = 7 
basins of the Mississippi. The linear correlation coeffl- 
cient is r = —0.06 and the Spearman correlation coeffi- 
cient is rg = —0.08. The latter has probability ps = 0.43 
indicating there are no significant correlations. 



C Effect of basin shape on Hack's law 

In what follows, we extract two statistical measures of 
correlations between deviations in Hack's law and overall 
basin shape. These are r, the standard linear correlation 
coefficient and r,,, the Spearman rank-order correlation 
p8|. For TV observations of data pairs 



coefficient 66, pi 
(ui,Vi), r is defined to be 



Y^f^iiui - fJ'u){vi - fJ.v) _C{u,v) 



(33) 



where C{u, v) is the covariance of the u^'s and w^'s, fiu and 
fly their means, and au and their standard deviations. 
The value of Spearman's is determined in the same 
way but for the Ui and Vi replaced by their ranks. From 
Tg, we determine a two-sided significance ps via Student's 
t-distribution Q. 

We define k, a measure of basin aspect ratio, as 



K = /a. 



(34) 



Long and narrow basins correspond to k ^ 1 while for 
short and wide basins, we have /t ^ 1 

We now examine the discrete derivatives of Hack's 
law in more detail. In order to discern correlations 



between large scale fluctuations within individual basins, 
we specifically look at the last two differences in a basin: 
^a,n-i and ft.n-i,n-2- For each of the Mississippi's 104 
order Vl = 7 basins, these values are plotted against each 
other in Figure |l^. Both our correlation measurements 
strongly suggest these differences are uncorrelated. The 
linear correlation coefficient is r = —0.06 ~ and, simi- 
larly, we have Ts = —0.08 ~ 0. The significance's — 0.43 
implies that the null hypothesis of uncorrelated data can- 
not be rejected. 

Thus, for Hack's law in an individual basin, large scale 
fluctuations are seen to be uncorrelated. However, corre- 
lations between these fluctuations and other factors may 
still exist. This leads us to our second test which con- 
cerns the relationship between trends in Hack's law and 
overall basin shape. 

Figure ^ shows a comparison of the aspect ratio k and 
/17.5 for the order Q. = 7 basins of the Mississippi. The 
measured correlation coefficients are r = 0.50 and = 
0.53, giving a significance of ps < 10~^. Furthermore, 
we find the differences hj^ (r = 0.34, = 0.39 and 
Ps < 10^^) and ft.6,5 {r — 0.35, Vg = 0.34 and Ps < 
10^^) are individually correlated with basin shape. We 
observe this correlation between basin shape and trends 
in Hack's law at large scales, namely ft.n,a_i, /ir2_i,r2-2 
and /if2,a-2, repeatedly in our other data sets. In some 
cases, correlations extend further to /ia_2 a-3- 



14 



0.65^ 

0.6 
0.55 

0.5 
0.45 

0.4 



-e- 



o 



o 



o 



o 



o 

o ° o o 

O - C r, 



O 



O 



() 




O 



0.2 0.4 0.6 0.^ 



log^Q(K = L /a) 



FIG. 13: Correlation between trends in Hack's law and 
the aspect ratio of a basin as estimated by k = L^/a. The 
data is for the order fl — 7 basins of the Mississippi and 
the specific trend is /ly 5. The correlation measurements 
give r = 0.50, = 0.53 and ps < 10~^. 



Since the area ratio Ra is typically in the range 4-5, 
Hack's law is affected by boundary conditions set by the 
geometry of the overall basin down to sub-basins one to 
two orders of magnitude smaller in area than the over- 
all basin. These deviations are present regardless of the 
absolute size of the overall basin. Furthermore, the origin 
of the basin boundaries being geologic or chance or both 



other basins. 

With the idealized Hack distribution proposed, we may 
begin to understand deviations from its form. As with 
any scaling law pertaining to a physical system, cutoffs in 
scaling must exist and need to be understood. For small 
scales, we have identified the presence of linear sub-basins 
as the source of an initial linear relation between area and 
stream length. At large scales, statistical fluctuations 
and geologic boundaries give rise to basins whose overall 
shape produces deviations in Hack's laws. Both devia- 
tions extend over a considerable range of areas as do the 
crossovers which link them to the region of intermediate 
scales, particularly the crossover from small scales. 

Finally, by focusing in detail on a few large-scale exam- 
ples networks, we have found evidence that river net- 
works do not belong to well defined universality classes. 
The relationship between basin area and stream length 
may be approximately, and in some cases very well, 
described by scaling laws but not exactly so. The gradual 
drift in exponents we observe suggests a more complicat- 
ed picture, one where subtle correlations between basin 
shape and geologic features are intrinsic to river network 
structure. 
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